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The extended t — J model is theoretically studied, in the context of hole underdoped cuprates. 
Based on results obtained by recent numerical studies, we identify the mean field state having 
both the antiferromagnetic and staggered flux resonating valence bond orders. The random-phase 
approximation is employed to analyze all the possible collective modes in this mean field state. 
In the static (Bardeen Cooper Schrieffer) limit justified in the weak coupling regime, we obtain 
the effective superconducting interaction between the doped holes at the small pockets located 
around k = (±7r/2, ±7r/2). In contrast to the spin-bag theory, which takes into acccount only the 
antiferromagnetic order, this effective force is pair breaking for the pairing without the nodes in 
each of the small hole pocket, and is canceled out to be very small for the d x 2_ y 2 pairing with nodes 
which is realized in the real cuprates. Therefore we conclude that no superconducting instability 
can occur when only the magnetic mechanism is considered. The relations of our work with other 
approaches are also discussed. 



I. INTRODUCTION 



Since the discovery of high temperature superconduc- 
ting cuprates^ it has been established that the strong 
Coulomb repulsion between electrons plays an essential 
role in the physics there^ and the magnetic mechanism 
of the superconductivity has been studied intensively. 
There are two streams of thoughts; one is the antifer- 
romagnetic (AF) spin fluctuation exchange based on the 
(nearly) AF ordered state, while the other is the reso- 
nating valence bond (RVB) mechanism with the focus 
being put on the spin singlet formation by the kinetic 
exchange interaction J. The representative of the for- 
mer is the spin-bag theory^ where the doped carriers 
into the spin density wave (SDW) state form small hole 
pockets, and exchange the AF spin fluctuation to result 
in the d xy superconductivity. The z-component of the 
spin fluctuation, i.e., \ zz \ gives the dominant contribu- 
tion. The RVB picture, on the other hand, puts more 
weight on the spin singlet formation and takes into ac- 
count the order parameter defined on the bond but usual- 
ly does not consider the antiferromagnetic long range or- 
dering (AFLRO )£ These two scenarios have been stu- 
died rather separately thus far, and the relation between 
them remains unclear. Hsu was the first to take into 
account both the AFLRO and the RVB correlation at 
half-filling in terms of the Gutzwillcr approximation. 6 
His picture is that the AFLRO occurs on top of the d- 
wave RVB or equivalently the flux state. It is also sup- 
ported by the variational Monte Carlo study on Heisen- 
berg model, which shows that the Gutzwiller-projected 
wavefunction ^rvb+sdw starting from the coexisting d- 
wave RVB and SDW mean field state gives an excellent 
agreement with the exact diagonalization concerning the 



ground state energy and staggered moment i The varia- 
tional wavefunction projecting the simple SDW state, 
on the other hand, gives higher energy. These results 
mean that both aspects, i.e., SDW and RVB, coexist in 
the Mott insulator. At finite hole doping concentration 
x, the AFLRO is rapidly suppressed and the supercon- 
ductivity emerges. Here an important question still re- 
mains, namely short range AF fluctuation dominates or 
the RVB correlation is more important. Because there is 
no AFLRO and no distinction between x zz and x ± , this 
question might appear an academic one with the reality 
being somewhere inbetween. Well-known results from the 
variational Monte Carlo studies are that the variational 
wavefunction V&rvb+sdw has the lowest energy therefore 
RVB (superconductivity) and SDW coexist also at small 
doping level£ From the viewpoint of the particle-hole 
SU(2) symmetry, this means that the degeneracy, i.e., 
the SU(2) gauge symmetry, between d-wave pairing and 
7r-flux state is lifted at finite x, and the former has lower 
energy. 

However this accepted view has been challenged by 
recent studies based on the exact diagonalization^ and 
variational Monte Carlo method^ These works found 
that the energetically most favorable state with small 
pockets at small hole doping levels is not superconduc- 
ting. This nonsupcrconducting state is consistent with 
the doped SDW+7r-flux state with the hole pockets near 
k = (n/2, ±7r/2). This means that the mean field picture 
is not so reliable at small x because the phase fluctuation 
of the superconductivity is huge at small x where the 
charge density, which is the canonical conjugate operator 
to the phase, is suppressed. Once the superconductivity 
is destroyed by this quantum fluctuation, the only order 
surviving is the AFLRO, and the system remains non- 
superconducting. Considering that the AFLRO state at 
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x = is well described by the SDW+(c?-wave RVB or 
7r-fhix) state, the relevant mean field state at small but 
finite x is SDW+7r-flux state with small hole pockets. 

From the experimental side, recent ARPES datas 
on Na-CCOC found the small "Fermi arc" near fc = 
(n/2, ±tt/2), and the fc-dependence of the "pseudo-gap" 
is quite different from that of (cosk x — cosk y ) expected 
for the d x 2_ y 2 pairingiifliii This result strongly suggests 
that the pseudo-gap is distinct from the superconduc- 
ting gap, and the superconductivity comes from some 
other interaction(s) different from J. Although only the 
"arc" is observed experimentally, the Fermi surface can 
not terminate at some fc-points inside of the Brillouin 
zone (BZ). Considering that the Fermi surface disap- 
pears with a large pseudogap at the anti-nodal direction, 
i.e., near fc = (jr, 0) and (0,tt), it is natural to assume 
that the small hole pocket is formed, along half of which 
the intensity is small and/or too broad to be observed 
experimentally^ Then the question of the pairing sym- 
metry arises because d x 2_ y 2 requires the nodes at the 
small hole pockets, which usually reduces the condensa- 
tion energy and is energetically unfavorable. The natural 
symmetry appears to be d xy without the nodes at the 
hole pockets, as has been claimed by the original spin- 
bag scenario^ 

In this paper, we study the superconductivity of the 
doped SDW+7r-fiux state for small x. This state includes 
both the RVB correlation and the AFLRO. There are two 
important effects of this RVB correlation; one is to intro- 
duce the parity anomaly to the nodal Dirac fcrmions at 
fc = (7r/2,±7r/2) and the other is to enhance the trans- 
verse spin-spin correlation x compared with the longi- 
tudinal \ zz . These two aspects might have crutial influ- 
ence on the superconductivity in the underdoped region. 
We have employed the 1 /TV-expansion, or random-phase 
approximation (RPA), to derive the effective interaction 
between the quasi-particles along the small hole pockets 
and the pairing force derived from it. 

The plan of this paper follows. In Section [H] we dis- 
cuss the model, formulation, and its mean field treat- 
ment. The Gaussian (second order) fluctuation around 
the mean field saddle point is treated and the effective 
interactions between the quasi-particles are studied in 
Section mil Section HVI is devoted to discussion and con- 
clusions. 



II. MODEL AND MEAN FIELD THEORY 
A. Formulation of the microscopic model 

The microscopic model of cuprates we consider for the 
study is the well-known t — J modeli^i^ The second 
and third nearest neighbor hopping terms are taken into 
account^ as they are necessary to describe the ARPES 
experiments measurements^ We use the slave bosons 



represcntationiSiii of the electronic operators 

4,. = flh, (1) 

with the constraint 

b\bi + E ftji,a = 1. (2) 

In terms of this formalism, the double occupancy of 
each site is excludedi 18 i 19 i 20 The operators f\ f are 
the fermionic ones while b\ b are bosonic ones (slave 
bosons) iSiS In all what follows we will exclusively work 
at zero temperature, and all the slave bosons will be as- 
sumed to be condensed. (In actual calculation, we take 
an extremely low temperature by technical reason. The 
results, however, are saturated and can be regarded as 
those at zero temperature.) 

We consider the following Hamiltonian on a square lat- 
tice with the lattice constant put to be unity and con- 
taining N s sites 

«= - (E*tf + £4-+ E 4) 

W> (i,0)' (i,3)" 

x E [ b i b ]flJj,°+ b AflJi,« 

cr=M 

+ ./E S '- S ^ ( 3 ) 

(i,3) 

where (), () , () denote nearest neighbor, second nea- 
rest neighbor and third nearest neighbor pair, respec- 
tively. We assume that the hopping contribution is uni- 
form: t i:j = t , t' tj = t , = t". For t' < , t" > 
the previous Hamiltonian deals with hole doped (p-type) 
cuprates, while the case t' > , t" < concerns elec- 
tron doped (n-type) cuprates. The unity of energy is: 
t = 0.4 cV = 1, and we will only study the hole doped 
case. Following Ref . Il5l we assume t = — 0.3, t — 0.2 for 
the extended hopping terms, and J = 0.3 for the Heisen- 
berg interaction. We have checked that slight variations 
in the numerical values of these parameters do not induce 
any significant modification of our results, as discussed 
in Section InTCl 

We construct the mean field theory to study this 
Hamiltonian. It has been recognized that the validity 
of the mean field theory is much less trivial in the theory 
with constraint. The simple comparison of energy does 
not offer the criterion because the mean field theory vio- 
lates the constraint, and hence can lower the energy in 
the physically forbidden Hilbert space. Therefore the 
choice of the mean field theory requires a physical in- 
tuition, and we employ the nonsuperconducting saddle 
point according to the reasons explained in the Intro- 
duction. In order to study simultaneously the influences 
of antiferromagnctism and flux in the resonating valence 
bond state, we introduce the two mean field parameters. 
The first one is a staggered magnetization^ 

(7,(7 
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FIG. 1: (Color online) Square lattice exhibiting simultane- 
ously a staggered magnetization and a flux phase. 



with s m = x, y, z, and a Sm being the Pauli matrices. 
The second one is the flux phase paramcte»2ii2i 

Xij = ( E flak") ' X = \Xijl 



Xi,i+x = X ■ exp 
Xi,i+y = X ■ exp 



(5) 



Fig. n shows the pattern of the flux and the staggered 
magnetization on the square lattice. 

We rewrite the model Hamiltonian J2J as 



with 



H = 



-E*-E*'- E*' 



(6) 



H m = aJ Y Si • Sj , 

= (l-a)J^S l -S J 



(7) 



The mean field Hamiltonian for each parameter is given 
by 



n 



MF 



<*J E E K'"Sf"+<"S 

s m =x,y,z 

(1 - a) J 



(1 - a) J 



[Xij fj, cr.fi, cr + Xij fi,cr fj.v] 

E Xij Xij- 

(id) 



In the previous expressions a parameter a has been in- 
troduced to divide the Hciscnbcrg exchange interaction 
into AF part and flux part. Here, the value of a is de- 
termined so as to reproduce the optimized result by the 
Gutzwiller variational method at half-filling. 6 Hence a 
should be regarded as a variational parameter, but we 
will keep it constant (a = 0.301127, as explained in Sec- 
tion lll C|) independently of the doping. We work in a path 
integral formalism by using the Stratonovitch - Hubbard 
transformation^^ 6 . The partition function is written 



V^V^VxDmVXDb exp 



dr £(r) , (8) 



E + bjb\fj i(T fi !Cr , with (3 = 1/ksT the thermal factor and the Lagrangian 



£ M = - (E*+ E *' + E E [«'i(r)6](r)* i , (T (r)* J> (r) + 6 J -(r)6j(r)* J> (r)* 4>(T (r) 
+ J2^iAr)[dr-Mr) - /i]* ilCT (r) +^&J(r)[5 T - iA^r)]^) 

+ « J EE E k s w*..(^:/,/M+<w^w?:,^vw] 

- E E [x«M*i,«rM*i,»(r) + xUr^iA^jAr)] 

(id) a 

- «J]T K(r)mJ(r) + m?(r) m ?(r) + m*(r)mf (r)] + ^i-^ £ Xy ( r ) X £(r). 



(9) 



In Eq. © "F^o-, are the Grassmann variables associated with the fj a , fi i(T operators, respectively. The Lagrange 
multipliers (iA^ ) assure that the constraint of no double occupancy is satisfied. The chemical potential \x associated 
with the fermions controls the electron density n/, and hence the hole dopping level x. 
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B. Saddle point hypothesis and derivation of the 
mean field Hamiltonian 



To identify the saddle point solution we consider the 
following assumptions. As we work at zero tempera- 
ture, all the holons are supposed to be condensed, i.e., 
{bi)o = bo with (bo) 2 = x where x is the hole doping 
concentration. Therefore the / operators can be sim- 
ply viewed as the renormalized electron operators. The 
Lagrange multipliers are considered independent of the 
sites, which signifies that the constraint is imposed on 
average at the mean field level 



(iAi)c 



An 



(10) 



For the antiferromagnetic part we impose an AF order 
polarized in z-direction 

(mf)o - , (mf)o - , (mf)o - (-l)*m. (11) 

Concerning the flux contribution we impose at half-filling 
(x = 0) a 7r-flux scheme (<f> = tt), which is energetically 



J 



the most favorable case for a square latticed 
{Xi,i+x)o = (Xi)o = X ■ exp 



(Xi,i+y)o = (xDo = X ■ exp [-i- (-1)' 1 



(12) 



To write the mean field Hamiltonian TImf, we note 
that the wavevector of the AF and 7r-flux is Q = (tt,tt). 
Therefore the summation over the momentum is done in 
the reduced magnetic first BZ. The borders of the ma- 
gnetic (or reduced) BZ are given by the Fermi surface 
of the Hubbard model at half filling. The reduced BZ is 
characterized by the nesting property under the vector 
Q = [tt,tt). TImf can be expressed in an appropriate 
spinor basis of the momentum space 



/fc,cr 
fk+Q,cr 



(13) 



We obtain the mean field t — J Hamiltonian written in a 
matrix form 



k+Q,a 



-AP -A* 



Cfc+Q 



fk,a 
fk+Q,< 



(14) 



r 



The summation ' in fc-spacc is over the reduced ma- 
gnetic BZ, and we have defined the following energies 

£fc = —2t[cos(k x ) + cos(ky)] — 4t [cos(k x ) . cos(ky)] 

-2t" [cos(2fc x ) + cos(2fc y )] , 
£fe = £fcX - (An + fi) 

-(1 - a) Jxcos(0/4) [cos(fc x ) + cos(k y )] , (15) 
Cfc+Q = £fc+QX - (A + n) 

+ (1 - a) Jxcos(0/4) [co8(k x ) + cos(k y )] . (16) 

The order parameters associated with the staggered ma- 
gnetization and </>-flux are respectively 

AJ^ = 2aJma , (17) 

A fc,<r = i(l - a)Jxsin(0/4)[cos(fc K ) - cos(fc y )] , 

for a = ±1. TImf can be diagonalized as 

kmf = £ ' E [KHk^i*,* + E l nka^\ , 



by a Bogoliubov-Valatin unitary transformation^^ 



fk,a 




Uk,a V k 




7lfe,cr 


fk+Q,a _ 








_ J2k,a _ 



(18) 



In the previous expression T lh<<r , 7ifc,o- (ltk,<T> 72fe,<7 ) 
are the creation and annihilation operators of the upper 



(lower) Hubbard band, respectively, with the correspon- 
ding energy eigenvalues (see Fig. EJ) 



K p 



ik + ik+Q 



(19) 



y[^-^ +Q ] 2 +A(\A^ + \Ai 

£,k + £k+Q 



*k,a\ 



(20) 

The elements of the unitary matrix are given by 

u k ,a = cos(9 k ) . e^" , v k = sin(6» fc ) , (21) 
with the trigonometric factors 
cos(6> fc ) 



\ 



1 



£fc - £fc+Q 



(a-a+Qr + 4(|A™J2 + | A x | 2) 



sin(6» fc ) 



\ 



(a-a+Qr+4(iA™j2 + | A x | 2) 
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&, (Femii level, hole doped) 



FIG. 2: (Color online) Upper and lower Hubbard bands asso- 
ciated with the eigenvalues E^ p and E]^ w of "Hmf- 
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FIG. 3: (Color online) Dependence of the modulus of the 
staggered magnetization m obtained at the saddle-point, as a 
function of the doping x = 1 — n/. 



Now wc expand the action up to the second order with 
respect to the deviation from the saddle point solution. 
At the first order we can derive the self-consistent mean 
field equations. The details are given in Appendix^ and 
we give here the final results 



sm(0 fc , 



A 



m 12 



!-^-"^y^' I uv 1 ; l=o. 



(24) 



It is noted here that the unitary transformation con- 
tains the complex phase factor e I ^ fc -' T , which becomes im- 
portant when we solve the BCS equation in Section ITTll 
In the field-theory terminology, this offers an example of 
"parity anomaly" in (2+l)Di2& 



-E 



[cos(fc a ) + cos(ky)] (£, k+Q - g fc )cos(0/4) 
2i[cos(fc x ) - cos(fcj,)]A£ CT sin(</>/4) 



(25) 



C. Saddle point solution and mean field equations 

In the path integral language, the saddle point action 



is 



So = ~ Yj 'EE [*fe,<r(iw„) *fc+Q.a(i^n)] 



x Q Q (fe,u,iw„) 

*fc,cr(i^ri) 



(22) 



with \io n the fermionic Matsubara frequencies and Qq the 
free Green's functions' matrix 

Go(k,a, iuj n ) 



-(iw„-a+Q) Aj^ + A*, 
Afc^ CT — A? -Qun-b) 



X 12 



(23) 



-E' 



[cos(k x ) + cos(ky)] (tk+Q - 6e) 



E u P 



= XCos(0/4) , 



r (Cfc+Q - £fc) 



E u P 



+(4 + t k ) 



(26) 



(27) 



= -An 



where the quantities 6h-Q> A* ^ E% p , E l k ow , t k , t k 
and t ' k are given in E qs. (T^l , (jPo|) . (|T7jl . ifT^I . J20J), 
(|Allj) . (|A12|) and (|A13j) . respectively. 

We adjust the electron number nj by the chemical po- 
tential (i. At zero temperature the upper Hubbard band 
is empty, while the lower band is partially filled as 



±y i 

N s ^ exp[f3E l ° w ] + 1 



1 -x, 



(28) 
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FIG. 4: (Color online) Dependence of the flux parameters (j> and x> Lagrange multiplier Ao and chemical potential fi obtained 
at the saddle point, as a function of the doping x = 1 — rif. 



where (J2k) 1S extended over the first BZ. 

We have solved numerically these equations by dis- 
cretizing the reduced BZ in 2 millions of points, providing 
us a precision on the obtained values better than 10~ 5 . 
In particular at half-filling (x = 0) we have found 

a = 0.301127, (29) 

by imposing the relation: m = 0.5% previously obtained 
by Hsupi Away from half-filling this relation between m 
and x wm be changed but the value of a <(2"9"|) is assumed 
to be the same on all the range of doping. This assump- 
tion means that the weight assigned to each of the de- 
coupling terms of the Heisenberg exchange interaction is 
kept constant as a function of the doping. This assump- 
tion does not change the essential features of our results 
presented below. 

We present the numerically obtained values of the 
mean field parameters as a function of the doping x in 
Figs. 01 and 0] We see that the Neel state disappears at 
a very small value of the doping, typically 1.5 %, which 
is in agreement with the experimental phase diagram of 
hole doped cuprates. This is in sharp contrast to the pre- 
vious studies of the t — J model, which found that the 
AF state could remain until doping values of around 15 
The flux 4> decreases relatively linearly from tt at 



half-filling, and reaches <\> ~ 2 for a doping of 10 %; this 
result is similar to the data obtained by Hsu et al.m^ 

We also show the shape of the Fermi surface for 
three doping cases: x = 0.01,0.05,0.1 in Fig. [SJ where 
the Fermi surface is composed of arcs which delimit 
the hole pockets located around the four nodes k = 
(±7r/2,±7r/2). 

In the next section we will consider the quantum 
fluctuation around this mean field solution. The second 
order Gaussian fluctuation corresponds to the RPA. We 
will calculate the "a la BC S"£2i21 pairing potential in 
terms of the exchange of these fluctuations, to see the 
effects of the coexistence of antifcrromagnctism and the 
staggered flux on the pairing force in undcrdoped region. 



III. GAUSSIAN FLUCTUATION AND BCS 
PAIRING INTERACTION BETWEEN 
QUASI-PARTICLES 

A. Second order action and correlation functions 

In this Section we expand the action with respect to the 
fluctuations of the mean field parameters up to second or- 
der starting from the saddle point solution for the t — J 
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model (|14p. By analogy with a diagrammatic language 
such a treatment is equivalent to taking into account 
all the bubbles associated with the fluctuating bosonic 
fields^ This approach allows us to calculate the diffe- 
rent correlation functions (or susceptibilities, both being 



linked via the fluctuation-dissipation theorem^) between 
those fields. 

Integrating over the fermionic Grassman variables and 
expanding the action with respect to the fluctuations of 
the bosonic modes up to the second order, we obtain 



S 2 = f drl - ]T(6 ) 2 [2 ■ iSXi ■ 6b, + A • (Sb,) 2 ] - aJ £ £ [5m?" ■ 5m s f] + J £ [6 Xij ■ 6 X *j] 

-(J2 t+ Y, t,+ 0(M'><^^-<%[(^*^} + (*^ . (30) 

(so) (i,i)" 

Evaluating the bubbles ^TrfCfoVi^oVi] as explained in Appendix 151 the quadratic action is given by 

S 2 = ^2 ' 5X l (q 1 ,iLj e ) (M i ,j(q,q 1 ,q 2 ,iwi) + (q, q l5 g 2 , iw £ ) j 5Xj(-q 2 , -iu; e ), (31) 

where the first order fluctuations of the bosonic fields SXi arc defined in Eq. 1|B5|I . and the matrix elements Mij and 
Hij are detailed in Eqs. (|B6|) and i|B12f) . respectively. In the previous formula, ioje denotes the bosonic Matsubara 
frequencies. 

In a path-integral formalism, the second order term of the effective action gives easy access to the different correlation 
functions between bosonic fields renormalizcd at a RPA level 35 



C(q,q 1 ,-q 2 ,we) 



l<i,j<9 



(Xi(q 1} iu)e)Xj(-q 2 , -iuji))rpa 



(X i (q 1 ,ioj£)X j (-q 2 , -VjJi))rpa 



l<ij"<9 



M(q,q 1 ,q 2 ,iuj e ) + -U(q,q 1 ,q 2 ,iujg] 



l<i,j<9 
1 



(32) 



The behaviour of the obtained transverse spin-spin correlation function \ ± is m agreement with well-known results 
concerning the Heisenberg antifcrromagnets, as discussed in Appendix ID1 
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B. Derivation of the effective action 

We can now explicitly calculate the "a la BCS " ??i? 4 
pairing potential exchanging all the collective modes 
described in the previous section. Following Schricffer 
et al.r^ we adapt their approach to the framework of 
the t — J model. We neglect the retardation effects spe- 
cific to the Eliashbcrg thcory22i2i and build a pairing po- 
tential by assuming the static limit (itu£ = 0). Therefore 
the frequency dependence will not be mentioned from 
now on. This is justified in the weak coupling region, 
where the adiabatic approximation kBT c ,Asc « f^D 
is satisfied ( T c : transition temperature, A$c'- supercon- 
ducting order parameter, u>d' frequency of the exchan- 
ging bosons). This approach is not sufficient to describe 
the superconducting state in underdoped cuprates where 
the features of the strong coupling effect are observed 
experimentally^ However, as will be shown below, the 
magnetic mechanism based on the generalized spin-bag 
theory gives only a very small pairing force or is pair 
breaking. Therefore this weak coupling approximation 
is justified a posteriori, even though it does not describe 
the real cuprates. 

We start with the linear interaction between the 
fermions and the bosonic fields (|A15|) 



*/ e =E'E 



^k+q,<j' ^fc+q+Q.er' 



x Vi(fe + q, a ; fe, a) 

^k+Q,a 



(33) 



Using the notations defined in Appendix iBl the interac- 
tion matrix Vi <|A16|) can be rewritten as 



Vi(fe + q, a ; fc, a) 



(34) 



E 



9 



1\=<1,<1+Q *i=l o'=l,2 



E </( fc '9i)V(9i) 



x«Q(c?i), 

where 8Xi, c i a > , s l a < arc given by Eqs. (jB5f) . {BSJ, (|B9fl . 
and IjBlOp . respectively. 

We remember that ty, ^ are the Grassmann variables 
associated with the f>, f spinon operators. As we are 
interested in the interactions between two 72 fermions of 
the lower Hubbard band, we introduce $, $ the Grass- 
mann variables associated with the 72, 72 operators. By 



neglecting the contribution of the upper Hubbard band 
we obtain from the diagonalization of the mean-field 
Hamiltonian 1)18(1 



* fcjCT = sin(6» fc )$ fci(T , 
*fe+Q,a = c-^ fc "cos(6l fc )$ fc! , 



(35) 



The first order action related to the 72 operators is given 
by 

St = E 'EW 

tj' ,(T 



sin(6» fc+ q) e C os(9 k+q ) 
9 



E E( E </(^i)*v(<*o 



q 1 .=q,q+Q i, = l a '=l,2 



x5Xi( qi ) 



x 



sm(U k ) 
e -^fc,-cos(6» fe ) 



(36) 



To incorporate the interaction effects at a RPA level, we 
build an effective action by adding to S( e the second 
order in <5JQ contribution £2 , i-e. using Eqs. (|31|l and 



s * = E ' E E 

1 Qi,i2=i^i+Q i .i= 1 



(37) 



-xC. l j(q,q 1 ,-q 2 )-SX 1 {-q 2 ), 
then the effective action S e ^ is given by 

S eff =S( e +S 2 . (38) 

In order to obtain the pairing potential, we integrate 
out the bosonic fields SXi in S e ^ . For convenience we 
define (^(q^ for q 1 = q, q + Q as 

<?(<7i) = [<M<7i)] 1<4 < 9 , 



(91) = ^E ' E [ sin (^+.) e i ^ 1 +-icos(e fel+g )] ( J2 <if 1(fel ' 9i)g ^i (9i) 

ki a 1 ,<ti a, —1.2 



shi(6l fcl ) 
e" i< ^ fc i.-icos(6' fci ; 



(39) 
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After integrating out the bosonic fields we have 



i qi,i2=<ii+Q 

or equivalcntly from Eq. (|39[l 

* e// = -jE' E EE' E 



xj sin(6» fc2 _q) e h 2—°2Cos(9 k2 - q ) 
x C i>j (q,q 1 , -q 2 ) 



sm(V k2 ) 

c -i<Afc 2 , CT2 cosf^fc 



sin((9 fcl+ q) e ^"icos^+q) 



E ^f^^v;^) 

0,1=1,2 



sin(6» fcl ) 
c-^ fc i-"icos(6' fcl ) 



(40) 



C. Calculation of the BCS pairing interaction 

To have a BCS form in the previous effective action 
(|4t)(l we impose the following relations 

* fcx = k , k 2 = — fc , q = k — fc, 

<r'i =T >o-' 2 =l, °i =1 , ax =T • (41) 

* fci = — fc . k 2 = k . q = k k , 

OX =1 , <72 =T , ^2 =T , (71 =| . (42) 

* fci = — k ,k 2 = k, q = k + k, 

a'x =T , ^ =| , (72 =T , <7i =| . (43) 

* fci = k , k 2 = —k , q = — fc fc, 

<r'i =i , °'<x =T , ^2 =4 , <7i =T • (44) 



The configurations (|41|l and (|42[) are related to the zz 
contribution, while the forms (|43[l and 144|) give the ± 
contribution. Therefore the BCS effective action is 

k,k' 

x *fc',T*-fc',i $ - fc ^ $ ^T- (45) 

where 

V BCS (k, fc') = VSb 8 {k, fc') + ^(fc, fc') (46) 

The zz BCS - type pairing potential V£q S is obtained by 
putting the configurations 1|41|) and l|42[) in the effective 
action <gU|) 



V^ cs (k,k') = £ £ 



q 1 ,q 2 =fe'-fc,fc'-fc+Q 1 

sin(0_ fe /) e I0 - fc '.icos(6 l _ fe ' 



X 



sin(0 fc 



a 9 — 1,2 

,.) e'V., cos(9fc ,)|( S c«( fc ,,0»„;( 5l ))[ e - 1 |:° ( lU)" 

0^=1,2 



sin(6» fc <) e 1 ^'.tcos(0 fe /) 



E c ^( fc .9a)2 i)a ;(9 a ) 

a 2 = l,2 



Sin(0 fe ) l r . /«._«.' __ n ) 
2- i ^T CO s(0 fc ) L ^ {K *' 9l ' ?2j 



sin(0_ fe 



) e i<#, -*'.icos(6L*/ 



i, =1,2 



(47) 
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FIG. 6: (Color online) Pairing potentials VgQ S {k,k ) 1471 1 . Vg CS (k, k ) 148H and VBCs(k,k ) 1461 1 as a function of k , with 
fixed k — (— 7r/2, — 7r/2), for two values of the doping: x = 0.01, 0.04. Numerical datas are plotted in unit of t. 



By imposing the configurations l|43[) and (|44|l in the effective action, we get the ± BCS - type pairing potential 



v± cs (k,k')= 

qi ,q 2 =k' +k,k' +k+Q M = l 



sin(0_ fe ,) e I0 -<='.icos(#_ fc O 
sin(0 fe /) c'^'.tcos^) 
sin(6*,') e^ h> icos(6. 



E c j;i< ^ '/.-,/•'/. 

i=l,2 



sin(0 fe ) 
e -i^ fc ,T cos(6> fc ) 

sin(0_ fc ) 
c" i ^- fc -icos(6'_fc) 



C itj (k +k,q 1 ,-q 2 ) 



a — 1,2 



sin(0_ fe ) 
e^- fc icos(6»_ fc ) 



sin(0_ fc O e I0 -'acos(0_ fc , 



sin(0 fc ) 
e -i<^,T C os(6» fc ) 



(48) 



The correlation functions and pairing interactions have fixed position k = (— tt/2, —tt/2), while the other one is 
been numerically calculated, the results of which are pre- at fc' moving inside the magnetic BZ. The intensities of 
sented in Figs. and These figures show the effective the effective interactions are plotted as a function of k : 
interactions between two 79 fermions, one located at the „ , , , , . . , ' , 

lor example the values given for k = {n/2, tt/2) corres- 
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FIG. 7: (Color online) Pairing potentials VgQ S {k,k ) 1471 1 . Vg CS (k, k ) 148H and VBCs(k,k ) 1461 1 as a function of k , with 
fixed k — (— 7r/2, — 7r/2), for two values of the doping: x = 0.06, 0.10. Numerical datas are plotted in unit of 



pond to a wavevector = fe — k = (tt,tt). As wg 
assume the static limit (iu>g = 0), the pairing poten- 
tials take real values. The dark color region means the 
repulsive interaction (V(k,k ) > 0), while the light re- 
gion displays the attractive interaction (V(k,k ) < 0). 
An interesting feature we can observe in Figs. and [7| 
is that the total pairing potential Vbcs is almost com- 
pletely driven by the zz (i.e. longitudinal or "spin-bag") 
contribution, without any major modication introduced 
by the ± (i.e. transverse or "RVB") correlation. As men- 
tioned in Section lll Al the results displayed in Figs. Eland 
[7| have been obtained by assuming t = — 0.3, t =0.2, 
and J = 0.3. We have also considered several additional 
sets of numerical values for the parameters in the range: 
-0.35 < t < -0.3, 0.2 < t" < 0.25, 0.3 < J < 0.35. We 
do not have observed any significant modification of the 
results. The robustness of these results with respect to 
slight variations of the parameter values can be explained 
as followed. In the present approach, the fermion spec- 
trum contains two Dirac fermions and associated small 
hole pockets at the nodal points. The appearance of 
Dirac fermions means that a topologically nontrivial fea- 
ture is introduced, which does not exist in the previous 
spin bag theory. Such a property gives qualitatively spe- 
cific characteristics to the system which are independent 



of slight variations of the parameter values, and therefore 
can explain the robustness of the results. In the next sec- 
tion, we will discuss the implications of the interaction 
Vbcs orL the possible superconductivity in this model. 



IV. DISCUSSION AND CONCLUSIONS 

In previous sections, we have calculated the effective 
interactions between the two holes in the background of 
the AF and RVB orders. There are four half-pockets of 
holes in the reduced first BZ as shown in Fig. 7, or the 
two hole pockets in the shifted BZ. These small hole po- 
ckets around k = (±ir/2, ±7r/2) are different from those 
of the simple SDW state as assumed in the original spin- 
bag theory^^ Namely the description of the antiferroma- 
gnetic state includes the singlet formation via the RVB 
order parameter in addition to the Neel order, and corres- 
pondingly the wavefunctions of the doped holes are dis- 
tinct from the SDW state. Already there opens the gap 
near k = (±7r,0), (0,±7r) due to the staggered flux, and 
the superconducting pairings near these points are irrele- 
vant. However the pairing force is still dominated by the 
zz-component of the spin susceptibility, and is repulsive 
for the momentum transfer = ±(7r, 7r). When the su- 
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perconducting gaps at k = (±7r, 0), (0, ±7r) are dominant, 
this leads to the d x 2_ y 2 pairing as discussed by several 
authora^i since the sign of the order parameter is oppo- 
site between fe = (±7r, 0) and k = (0,±7r). However in 
the present case, only the states near k = (±ir/2, ±7r/2) 
could contribute to the pairing. One of the main predic- 
tions of the spin-bag theory&a is that the superconductiv- 
ity occurs through the pairing of holes in the small pock- 
ets, which are delimited by the Fermi surface as we can 
see in Fig.|SJ The present study shows that the total pair- 
ing potential Vbgs is clearly repulsive, or pair breaking, 
on all the range of doping for a wavevector qd in the (0, 0), 
(ir, 0) and (0, ir) areas of the magnetic BZ. An attractive 
behaviour is also found when the momentum transfer is 
around (tt,tt). It corresponds to relatively large values of 
qd- Namely, the repulsive interaction (dark contribution) 
near k = (— ir/2,~- ir/2) is dominating over the attractive 
(light) region with large qd- Therefore as long as one 
considers the pairing form which is constant over each of 
the small hole pocket, there occurs no superconducting 
instability because the intra-pocket pair breaking force 
is larger than the inter-pocket force independently of the 
relative sign of the order parameters between different 
pockets. 

On the other hand, one can consider the pairing with 
the nodes in the small hole pockets as in the case of 
d x 2_ y 2 pairing realized in real cuprates^2ii£ In this case, 
however, most of the interaction cancels within the small 
pocket due to the sign change of the order parameter. 
We could not determine the sign of the residual pairing 
interaction, but we can safely conclude that it is very 
weak even though pair creating. Therefore we conclude 
that the doped staggered flux state is stable against 
the superconducting instability when only the magnetic 
mechanism is considered. This is in accordance with the 
recent exact diagonalization studies on the t— J model at 
small doping»£ Furthermore this suggests that the other 
mechanism of superconductivity is active in cuprates at 
least for small hole doping region. Also a related work has 
been done by Singh and Tesanovic^ where the quantum 
correction to the spin-bag mean field theory is studied 
up to one loop level and also the doped case is consid- 
ered. This is along the same line as the present study, 
and they also obtained the repulsive interaction between 
two holes. The new aspect introduced in our study is the 
RVB correlation represented by the flux order parameter. 
Combining these two works, the pair breaking nature of 
the magnetic interactions seems to be rather robust in 
the low hole doping limit. 

There are two possible routes leading to the super- 
conductivity. One possibility is that the starting mean 
field ansatz becomes inappropriate for optimal doping 
case. According to the SU(2) formulation of the RVB 
state, the staggered flux state is the quantum mechani- 
cal mixture of the flux state we discussed above and the 
d x 2_ y 2 pairing state, and the latter one is more and more 
weighted as the doping proceeds. Therefore our formu- 



lation is valid only for the small doping region and can 
not capture the crossover to the superconducting state 
since it is based on the perturbative method around the 
AF + flux state. On the other hand, starting from the d- 
wave superconducting state in the SU(2) formalism, the 
staggered flux fluctuation has been studied42il4 This is 
another language describing the instability towards the 
AF ordering, which is represented by the chiral symmetry 
breaking in the gauge theory. It is also found the self- 
energy due to the staggered flux fluctuation leads to the 
quasi-particle damping strongly anisotropic in the mo- 
mentum spaced These works are the approach from the 
superconducting side, and are complementary to ours. 

Taking the view that the magnetic interaction is pair 
breaking in the low hole doping limit, we should look 
for other forces which are active in the cuprates. One 
of the most promising candidates for this is the electron- 
phonon interaction^*^ which already manifests itself in 
the angle-resolved photoemission spectroscopy^ Howe- 
ver, much more work is needed to establish this scenario. 
Especially the interplay between the strong correlation 
and the electron-phonon interaction remains an impor- 
tant issue to be studied. 

In conclusion, we have studied the extended spin-bag 
scenario taking into account the staggered flux resonating 
valence bond order in addition to the antifcrromagnctic 
order. The pairing potential between the quasi-particles 
has been calculated at the Gaussian level, and it is found 
that the magnetically mediated interaction is pair brea- 
king or very weak at the small hole doping limit, sug- 
gesting the other mechanisms such as electron-phonon 
interaction arc active. 
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APPENDIX A: DERIVATION OF THE 
SELF-CONSISTENT EQUATIONS 

We allow the different fields to fluctuate at the first 
order around their saddle point values l|lUfl , i|ll|) an( l 11211 

bi = b (l + Sbi) ; 
iA ; = A + iSX, ; 

to-" 1 = (TO- m ) + <5m? m , s m = x,y, z ; 

xm+ Sx = xt x = (xt x )o + $xt x , 

Sxt x = S Xi x + ^ Xi* > s x = x, y. 
By expanding the fluctuations at the first order we obtain 
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+ j E [(*«)° • + (*«)s ■ fan] - E i iSX * ■ *v*<,*] 
- (E i+ E*' + E *")E( 6 o) 3 (<»i+^)[*i,^+* J >*i,a] 

W> («)" 

+ f J EE E [K m (^^:>^.o+^-(* 3> ^* ji(r o] 

(»,j> <t,<t' s m =x,y,z 



In the preceding formula, like in those which will follow 
in the appendices, the imaginary time dependence of the 
Grassmann variables and first order bosonic fluctuations 
follows the Lagrangian © and is implicit. For clarity 
in the following calculations we decompose S\ into two 
parts 



Si = 5f° + S( 



fe 



(A2) 



Si contains terms with exclusively auxiliary bosonic 
fields and no Grassmann variables. It corresponds to the 
three first terms of Eq. (|A1|) 

Si 10 = drj - Y^(bof(iSXi + 2A • 5h) (A3) 
—aJ J2 [(-1)'™ • <5mj + (-l)J'm • 8mf] 

+ 1 9 ~ J E [fei)o • $xtj + (Xij)o ■ s Xij] 



S[ e includes all the terms containing fermionic variables 
It takes the four last terms of Eq. (|AT|) 



S[ e = [ dri - V [i<5A, • 

-(E £ +E £ '+EO(M 2 

<*j>" 



(A4) 



+f J EE E [K"(*<^>«v) 

(t,3> o-.cr' s m =x,y,x 

We start by focusing on the bosonic part iS^ of the 
action Si expanded at the first order in fluctuations. By 
passing in the space of moments and Matsubara frequen- 
cies we obtain 



5{° = ~ (bo) 2 iS\ q=0 (iLu e = 0) (A5) 

-2(fe ) 2 A ■ 5b q=0 {iLu e = 0) 
+4aJ ■ m ■ 6m q —Q(iu>£ = 0)] 

+(i-a)jJ2 [x v - q {-m)-6x5(m) 



+X-q(-i^) • S X y q {iu e )], 

where iui labels the bosonic Matsubara frequencies and 
is expanded over the first BZ. 

We study now the fermionic contribution S[ e and de- 
compose it into four contributions 

S( e = S? + 5f + S[ n + 5f . ( A6) 

By passing in the space of moments and Matsubara fre- 
quencies we obtain 
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E ' 



- lUJ m ,lUJ n 



-i(5Aq(iw m - iw„) -i(5A q+ Q(iw m - iw n ) 
4SX q+Q (iuj m - iuj n ) ~\5X q (iuj m - iw n ) 



(A7) 



*fe,cr(i^n) 
^k+Q,a{^n) 



5 f = ^p=E 'E E 4,4W( iw ™)^+ g+ Qy( iw ™) 



\2 



(A8) 



-2(6 r(ife + ife+q + t k + t k+q + t k + t k+q )5b q (iu m - iu n ) 

-2(6o) 2 (-t\+ t k +q + h + t k+q + i' k + i' k+q )5b q+Q (i 
~2(b ) 2 (t k - t k +q + i'k + i'k+ q + i'k + tk+ q ) Sb q+Q( iuj ™ ~ i^n) 



lW m - VjJ n I 



-2(&0) (~tk - tk+q +t k + tk+q + l k + t fc+q) <5& q( iw " l ~ iw„) 



*fe+Q,<r(iW n ) 



E 'E E E ^Jv^'^J^+^+o,^ 1 



k,q 



aJ q • <Jm* m (iu; m - iu n ) -aJ q • £m*™Q(iu;, 
-aJ q ■ 5m^g(iu; m - kj n ) aJ q • (5m*™ (ia; 



771 1^71 J 



(A9) 



*fc,cr(iw„) 

*fe+Q,<r(iw n ) 



^ = TPrE'E E E °l',c [*k+ q ,*'(}"m) *k+ q +Q,*'(Mm) 

■Q.J2) 



(A10) 



(1 -a)J[ 



-S 



5 Xq (i^m - i^„) • COs(fe Sx , y Sx/ 

Xq x (iw m -iw n ) ■ sin(A: Sx + g^/2)] 



(1 - a) J[ - S'x q x +Q (iu m - iw n ) ■ sin(fc 

+ 6 "x q + Q {^ 

(1 - a) J [ + 6 Xq+Q^n - ■ sin(fc Sx + g Sx /2) 



Xq+Q^m - i^n) ' COs(fc 

' /ft f t ; 

x W 2 ) 
> m - iuj n ) ■ cos(fc Sx + q Sx /2)] 

r 



■9- x /2) 
-&*/2)] 



(1 - a)j[ + S'x q x (iuj m - ioj n ) ■ cos(k Sx + q Sx /2) 
""x q x (iu m - ioj n ) ■ sm(k Sx + q Sx /2)] 



*fc : <r(iw„) 

*fe+Q,<r(iw„) 



In the previous formulas cr°, is an element of the (2 x 2) 
identity matrix, iuj m is (like iw„) a fermionic Matsubara 
frequency and we have defined the following quantities 

t k = t[cos(k x ) + cos(ky)] , (All) 
t h = t [2.cos(fc x ).cos(fc y )] , (A12) 



t 



- t [cos(2fc x ) +cos(2fc y )] , (A13) 
J q = j[cos(q x ) + cos(q y )] . (A14) 



By using the Pauli matrices, we can put the different con- 
tributions HA7|> - (|A10|) together, and rewrite S{ e (|A"6|) 
as 



^ = E'E E 



[* k+q,a' ( icJ m) ^t+g+Q/K) Vi (fc + q, cr' , iw m ; fc, cr, iw„) 



*fc,<r(iw„) 

*fc+Q,<r(iw„) 



(A15) 
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with Vi the interaction matrix 



V^Ns Vi(fe + g,<r ,kj m ; k,a,kj n ) (A16) 
= E 1 [-(l-a)J ■cos(k Sx +q Sx /2)-a° (r , iT -a z ]S'x s q x (iuj m -iuj n ) 

s x =x,y 

+ [ - (1 - a) J ■ sm(k Sx + q Sx /2) ■ <j°, a ■ ia y ] 6 x^+qO^m - iw„) 
+ [ - (1 - a) J ■ sm(k Sx + q Sx /2) ■ cr°, a ■ a z ] 6 \q x - i^n) 

+ [+ (1 - a) J ■ cos(fc Sx +<? Sx /2) • tr°, )0 . • ia y ]S" x s q x +Q (iw m - iw„)| 

s m —x,y,z ^ ' 

+ | [ - <r°, a ■ a°]iSX q (iu} m - iw n ) + [ — o"°' jC . ■ o- :r ]i(5Aq + Q(iw m - icj n )| 

+ | [ - 2(6 ) 2 (i fe + + ti' + Cf<>°V ' d " ~ 2 ( fo o) 2 (4 + *fe+>° v • 5*]6b q (kj m - kj n ) 

+ [- 2(M 2 (i fc + i k+q + + C+gKv ' 5* + 2(6o) 2 (4 - 4+<>° v • iaf]5& q+q (iw TO - iw n )j . 

Finally the mean field equations are obtained after having integrated out the fermionic fields by checking the saddle 
point property^ 

Tr^oVi] = 0, (A17) 
which gives by using the previously obtained expressions J22J, l|A2j) . (|A5() . ()A15|) the set of coupled mean field equations 

APPENDIX B: DERIVATION OF THE GAUSSIAN ACTION 

We decompose the second order effective action 52 <|30|) into two parts 

S 2 =S b 2 °+S b 2 ub . (Bl) 

S 2 ° contains only contributions due to auxiliary bosonic fields. It includes all the terms of Eq. 1)30(1 with the exception 
of the trace 



S b 2 ° = f drl 5> ) 2 [2 • iSK ■ 8h + A • (Sb t ) 2 } - a J £ £ [Sm°™ ■ 5m s f] + ^^J^ [5 X « ■ Sxtj] 

J ° ^ i s m =x,y,z 

-( E 1 + E *' + E E(M 2< ^ • ^ [<*«.** 3 >) + <?i,**i,*)] } 



s m =x,y,z (ij) 

(B2) 



(<>J> <i,i>' <ij)' 



takes into account the "bubble" contributions 



S^ = ±Tr[g V l g V 1 ]. (B3) 
We firstly evaluate iS>2° <|B2(1 by passing in the space of moments and Matsubara frequencies 
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S h 2° = -(b ) 2 E ' E { 2i [ sx i( iu e) ■ Sb- q (-m) + 5\ q+q (iui e ) ■ 5b- q - Q (-m)] 



q 



+A [ Sb q (iuje) ■ 6b- q (—iwi) + 8b q+Q (iujt) ■ <$6_ g _Q(— iu^)] 
-«E ' E ^ [ 5m* m (io^) • Sm^i-iue) - 5m q l Q (iu e ) ■ 5m s ™ q _ Q {-\w t )] 



q vjJi s m —x,y,z 



q iui £ s x =x,y 



+ 5 'Xq+q(i^) ' ^xl x g_ q (-i^) + <S"x£_ Q (iw*) • J' x'* q - Q (-m)] 

t ■ l'f [cos(fe) + cos(q a )] + 2f ■ l%' v cos(q x ) ■ cos(g y ) + t ■ x [cos(2^) + cos(2g. y )] 

xSb q (iuj e ) ■ 5b_ q (-\LUi) 
+ < — t ■ l'i[cos(q x ) + cos(q y )] + 2t ■ l%' v cos(q x ) ■ cos(q y ) + t ■ 1% ' x [cos(2g x ) + cos(2q y )] 



x8b q+Q {\uji) ■ 5b_ q ^ Q (-iuj e ) 



where J q is given by Eq. (|A14|I and we have defined for s t , s t > = 



(B4) 



x,y 



i+ St +s t , ,a) — ^ 2-^i gup _ E iow cos \ K s t T «s 



In order to write S^ b (|B3|l and 6>2° (|B4fl in a compact way which will be useful in the following calculations we define 
the vector of nine components SX containing the first order fluctuations of the bosonic fields 



SX(q,iuj e ) = [5Xi(q,hJe)] 1<i<g = 



S''x q (^£) 
S"x v q {iuji) 
8m x q (\u)i) 
&m y q (\u)t) 
8m z q (iu)t) 

5b q (\u>t) 



Using Eq. I|B4|I we define the matrix M so that 



5X(q + Q, \u t ) = [SXi(q + Q, lut)] 1<i<g = 



5 X q+ Qfae) 

5 'Xq+Q(m) 

5 " n x q+Q {^i) 

S"x v q+Q ^l) 
8m x +Q (iujt) 
8m? q+Q (iuje) 

5\ q+Q {iuji) 
5b q+Q (\u)i) 



(B5) 



r 



•^(9,91,92)^) 



M i j(q,q 1 ,q 2 ,iwi) 



l<ij"<9 



5 2° = E' E EE^i'^) 

xMij(q, q 1 ,q 2 ,idj e ) ■ SXj(-q 2 ,-iwt) ■ (B6) 



We now tackle the calculation of the "bubble" part 
(lB3l). We haveS 



Tr^oV^oVi] (B7) 
E E E tr [ Go(k + q, a ,iw„ 

k q a, a' iw n,i^ 

xVi(k + q,a , iu n + iut, k, a, iw n ) 
x g (k,a,ko n ) 

xVi(fe, a,iui n ;k + q, a ,iw n + iui)] , 
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where tr is related to the spinor basis l|13(l . We define the 

vector c a ,a containing the coefficients associated with 
the first order fluctuations in the V\ matrix. As we can 
see in Eq. I|A16|) Sb q (and 5b q+ q) gives two separated con- 
tributions and not one like all the other auxiliary bosonic 



fields whose first order fluctuations are contained in the 
vector 5X <|B5[) . Therefore c a ,<J has ten components and 
not nine, which is indicated by the index a being equal 
to one for all the fields except 5b q (and 5b q +Q) for which 
it takes the values one and two 



<i<9 



-(1 

-(1 
-(1 
-(1 



a)Jcos(k x + q x /2)a°, 
a)Jcos(k y + q y /2)a°, 
a)Jsm(k x +q x /2)<j° cr ,[ 
a)Jsm(ky+q y /2)a°, 



a,J a a x , 
aJ a <7 v , 
aJ a a z , 



(B8) 



-2(bo) 2 (t k + t k+q +t k +t k+q )* a , 
-2(b ) 2 (t k +t k+q )a°, 



c° '°(k 7 q + Q) = [clf(k 7 q + Q)] 1 



<i<9 



- a)Jsm{k x + q x /2)a° a , ^ 
-i(l - a)Jsin(k y + q y /2)a°,' 
i(l-a)Jco8(k x + q x /2)a a ,' o 
i(l - a)Jcos(k y + q y /2)a° a , o 
— aJ a a x , 
—aJ a a y , 
—aJ„a z , 

H a ,<r 



-2(6 ) 2 (i fc + i k+q + t k + t k+q )*° a , a 
2i(b ) 2 (t k -t k+q )<T a , 



(B9) 



where the terms ffc, t k7 i k and J q arc defined in Eqs. i|All|) . I|A12|) . (| A13|) and l|A14(l . respectively. In the same 
manner we define the vector s to associate each fluctuation clement SXi to its corresponding Pauli matrix appearing 
in the expression of the Vi matrix (|A1B|> 



s(q) = [V («)] i< 4<9 



(BIO) 



The summation over the fermionic Matsubara frequencies iuj n is incorporated in 

F l ^(q 1 ),s ja „(q 2 )( fc '9,i^) = ^o(fe + g J ^\iw„+i^)s j)0 '(g 1 )^o(fe,o-,iw„)s 3 - »(g 2 ) ] , (Bll) 



with q 1 , q 2 = q, q + Q. Considering the identity matrix and the three SU{2) Pauli matrices we get from Eq. (|B11|I 
sixteen different thermal integrals which arc detailed in Appendix IU1 



Using the preceding expressions i|B5|l . i|B8(l . (|B9fl . 

(IBlOfl . and the formulas associated with the F - integrals, 
we define the matrix II which explicitly incorporates at 
a RPA ( "bubble" ) level the contributions of the different 
fields fluctuations to the correlation functions 



n(<?,gi,<z 2 ,i^) 



(B12) 



l<i,j<9 



with the matricial elements 



1 



E'E E </(Mi) 

« cr ,(T a ,a —1 

x ^;T(, l)tSj>o „ fa a ) (fc ' 9 ' ^ )c i,a" (fc ' 92) • 

From Eqs. l|B7)l. (|BTT|l and <fB12f> we get 

Tr[0 o ViGoVi] (BI3) 

= E' E EE 

<? qi,q2=q>q+Q ^ *j =1 

xll i!j (q,q 1 ,q 2 ,iui)5X j (-q 2 ,-iui) . 



With Eqs. (JBTJ, IJB6J, (B3J) and ljBT3|) we finally obtain 
<S 2 given by Eq. (|3*TJ) . 



APPENDIX C: THERMAL INTEGRALS 
EXPRESSIONS 

We give here the analytic expressions of the sixteen 
different F integrals (|B1 l|l useful to calculate at a RPA 
level the two (bosonic) fields correlation functions 1)32( 1 . 
We define 



Ki = [K P +q K. + i^)] [(iw„ + m) E l k ™ q ] 

x[E^-kj n ][kj n -Eir] , 

,9o,o = 2(iw„) 2 - io^(£ fe + Cfc+Q) 

+iw„ [2io^ - (£fe + £fc+Q + £fc+q + Cfe+q+Q)] , 

.go,. = -2io,„(A™ CT + A- w ) - 2i^A™ CT 

+ ^fc!cr(Cfc+q + ^fe+q+Q) + ^fcVq.o-' ^ fc + 6c+q) , 

9o, y = -2iu; B (iA* „ + iAj^y) - 2iwMJ i<r 

50,z = iw„(£fc — Cfe+Q + £fc+q _ ^fc+q+Q) 

(£fe - Cfc+Q) - (CfeCfc+q - Cfc+QCfc+q+Q) i 



— i^fc + i£fc+Q + i^fc+q - i^k+q+Q) 
— i^f(iCfe — i^fc+Q) + (i£fc£fc+q+Q — i^fc+QCfc+q) 

fez = 2iu; B (A* a - Al +g a .) + 2k* A*, 

~ A fc l0 -(Cfe+q + Cfc+q+Q) + ^k+q.cr' ^ k + ^ fc +Q) 

ff,,, = 2ia; n (iA^ - iA™ q></ ) + 2i^iA™ CT 

- iA fc"<j(Cfe+q + Cfc+q+Q) + lAfc+q,,/ (6c + Cfc+Q 

and we obtain 

F~o^o (k, q, iut) = ^ E i?3~ 1 5o,o 



+ (£fc£fc+q + ^k+QCk+q+Q) 

+2(AT AT , - A* A x ,) I 



F^(k,q,iuj e ) = - ^ * J g °° 



+ (£fc£fc+q+Q + Cfc+Q^fc+q) 

+2(A^A»VqV+A^A^,)} 



P J,., ^uA 



+ {£k£k+q+Q + ^k+Q^k+q) 

-2(AT AT , + Al A x 

\ fe,<7 fc+q,o- ' k,a k+q,cr 



F? m % (k, q, iu t ) = -3 ^ -=j- < 

+ (£fc£fc+q + ^fe+QCfe+q+Q) 

-2(A^A»VqV-A^A x w )}, 
^,5- 9' = 4 E pa- 1 go ' x 

"Afc^^fc+q - ^fc+q+Q) + A* , (£fe - £,k+ 



Q) 



q,iw t ) = ^E^ 



/3 



^o;^ y (fe, 9, m) = \ J2 pJ~ { so,,/ 

-iAfc'crlCfc+q - Cfe+q+Q) + ' l ^k +q , a ' (Cfe ~ Cfe+Q)| 
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%. (fe, 9, iw/) = ^ E £^ 1 ^ 



-2(AtT A x , - A? At™ ,) 



90, 



-2(AT A x , - A x At™ 
F^ a ay (k,q,iuj e ) = i^2-^r\9x,y 



+ 2i (A^A x +gV+ A x (T A™ +qy ) 



+2i(A^A x w + A x ^A™^,) 
F^(k,q,iuj e ) = I^-LL^ 



-y— 



^i 1 



-^<r(Cfc+q - Cfc+q+Q) - A fc+q lC r' (£fc ~ £fc+C?) 
- iA fc !( r(Cfe + q - Cfc+q + Q) - iA fc + (? y (& - Cfc+Q)| , 

^.;?v(fe,9,^) = 4E 

- iA fc, CT (^fe + q -^+q + Q) - iA fc + (? , CT ' -^+Q)} • 



APPENDIX D: EXISTENCE OF A GOLDSTONE 
MODE 



In this appendix we examine the contribution to the 
second order action Si given by rriq + Q and m^ + g in the 
half filled, uniform and static limit 



x = 0, q = 0, icj£ = 0. 



We write it as a part of 6>2 131|) 



(Dl) 



°2 



\^M 5 , 5 (q = 0, Q, Q, iut = 0) + in 5 ,5(9 = 0, Q, Q, \u t = 0) j [6m%(iu>t = 0)] 
\M 6 , 6 (q = 0, Q, Q, iiu e = 0) + in 6i6 (q = 0, Q, Q, iw* = 0)^J [5m y Q (h) t = 0)] 



« Jq=0 (<^TOq) + a Jg =0 (^q) 



^- E ' E [ - nJ " »<'•.„/" - [ - aJ <j=°< CT '] OK?) 2 - 



(D2) 



I 

By using Eqs. (|B6(I and (|B12|) we get in the simplified The mean field equation (|24|l related to the magnctiza- 
limit (|D1|) tion gives 

On 



2a J 



E 



, 8a 2 J 2 



-E 



, 8a 2 J 2 



N s ^ - 

fc fe « 



2aJ, 



(^™q) 2 + (<$™q) 



(D3) 



then with Eq. (|D3|) we have finally: 52 
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The behaviour of the second order action concerning 
the x and y directions of the magnetization implies that 
the transverse spin-spin correlation function x ± contains 
a gapless pole. It is a consequence of the Goldstone 
theorem^ which has to be applied in the present case be- 
cause of the rotational symmetry breaking in spin space 



due to the imposed AF order. It is in agreement with 
the effective field theory of quantum antiferromagnets, 
which was built by Haldane for one-dimensional chains^ 
and extended later to square lattice systems^ The exis- 
tence of this Goldstone mode has already been observed 
in the original spin-bag approach^*^ 



* Present adress: Max-Planck-Institut fur Festkorperfor- 
schung, Heisenbergstrasse 1, D- 70569 Stuttgart, Germany. 
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